%{
1.控制系统稳定的充要条件：闭环传递函数的极点全部在S平面的左半平面，即均具有负实部
2.判定稳定系统方法：
    1.直接判定法
    2.劳斯判定法：
        1）特征方程的全部洗漱为正值
        2）其劳斯表的第一列系数为正值
        特殊情况：某行第一个元素为0
        处理方法： 用小正数代替，计算劳斯表，不变号则系统临界稳定    
    3.波特图
        margin()--1）幅频响应：Pm:增益曲线下降1db时相位变化值
                  2）相频响应：Gm:相位变化值为180时增益与0db的距离
    4.奈奎斯特稳定性判断：
        1.. 首先确认开环传递函数右半平面极点个数P，可以直接通过开环传递函数得出
        2.. 绘制出开环传递函数 T(s) 的Nyquist 曲线（幅相特性曲线），查
            看其绕（-1，0）点的圈数 R，顺时针
            绕为负，逆时针绕为正
        3.. 通过 P-Z=R 得出 Z 的值，如果 Z=0
            则闭环系统不存在右半平面极点，系
            统稳定。反之系统不稳定
        奈奎斯特曲线绘制：nyquist(G);
%}
%%劳斯表创建
p=input('p');%系数矩阵
n=size(p,2);
if mod(n,2)==0
    n1=n/2;
end
if mod(n,2)==1
    n1=(n+1)/2;
    p=[p,0];
end
l=reshape(p,2,n1);
L=zeros(n,n1);
L(1:2,:)=l;
for i=1:2
       L(i,1)=max(L(i,1),0.001); 
end
for i=3:n
    for j=1:n1-1
        L(i,j)=(L(i-1,j)*L(i-2,j+1)-L(i-1,j+1)*L(i-2,j))/L(i-1,j);
    end
end
%%判断是否稳定
flag=1;
for i=1:n
    for j=1:n1
        if L(i,j)<0
            flag=0;
        end
    end
end
if flag==1
    fprintf("稳定");
else
    fprintf("不稳定");
end